Mode transformations and entanglement relativity in bipartite gaussian states 
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A proper choice of subsystems for a system of identical particles e.g., bosons, is provided by second-quantized 
modes i.e., creation/annihilation operators. Here we investigate how the entanglement properties of bipartite 
gaussian states of bosons change when modes are changed by means of unitary, number conserving, Bogo- 
lioubov transformations. This set of "virtual" bi-partitions is then finite-dimensionally parametrized and one 
can quantitatively address relevant questions such as the determination of the minimal and maximal available 
entanglement. In particular, we show that in the class of bipartite gaussian states there are states which remain 
separable for every possible modes redefinition, while do not exist states which remain entangled for every 
possible modes redefinition 



I. INTRODUCTION 

The physical notion of quantum entanglement lies at the 
very heart of the novel and rapidly developing discipline of 
quantum information science The mathematical defini- 
tion of entanglement requires the existence of a tensor product 
structure over the Hilbert state-space of the system S under in- 
vestigation; physically a multi-partition of S into subsystems 
has to be identified. It is nearly obvious that this latter choice 
is by far not unique; in most of the cases there is actually a 
enormous set of potential multi-partitions available. Unless 
a physical, possibly problem dependent selection criterion is 
somehow provided, all those "virtual" multi-partitions of S 
into subsystems are conceptually on the same foot. 

In Refs 121 and yf] an operationally motivated framework 
for relativity of entanglement, based on operator-algebras of 
observables has been introduced. There, the key physical 
idea is that the relevant multi-partition is selected on the ba- 
sis of the set of implementable quantum dynamics, includ- 
ing both unitaries and measurements, in the actual operational 
scenario. Mathematically, the necessary tensor product struc- 
ture is induced over the global set space by a bunch of imple- 
mentable commuting algebras of observables satisfying suit- 
able natural properties 1 3 ] . Each of this algebras becomes then 
a local one i.e., it plays the role of the observable algebras of 
one of the "virtual" subsystems. 

A very important class of physical systems is provided by 
the ensembles of identical particles. For this kind of systems, 
in view of the (anti)symmetrization postulates used to enforce 
the proper statistical behaviour, a tensor product structure of 
the physical state-space e.g., the totally symmetric subspace 
for bosons, is not naturally given. Accordingly, the extension 
of the notion of quantum entanglement to systems of indistin- 
guishable particle is not a totally trivial task and it indeed gave 
raise to different, complimentary approaches H, H. In this 
paper we shall follow the second-quantized approach where 
the subsystems are associated to modes |4], rather than to par- 



* ciancio@isiosf.isi.it 



1 giorda@isiosf.isi.it 

~~MJtl 



-! zanardi@isiosf.isi.i 



tides 0]. We will focus on bi-partite bosonic gaussian states 
and we will study how entanglement changes under unitary, 
number-conserving, redefinition of modes i.e., subsystems. 
Gaussian states play a central role in the theory of entangle- 
ment for continuous variables systems; the reason is two-fold. 
First their experimental relevance; indeed most of the states of 
light, which are easily obtainable in the lab by means of linear 
optics devices e.g., beam splitters, are gaussian. Second, they 
are mathematically easy to handle. In spite of the fact that 
gaussian states live in the infinite-dimensional bosonic Fock 
space all of their physical features are indeed captured by the 
so-called Covariance Matrix V. From this finite-dimensional 
matrix one is able to extract straightforwardly various mea- 
sures of entanglement like negativity and logarithmic negativ- 
ity 1 6]. To this aim it is often useful to cast the covariance 
matrix in standard form V s ■ This is possible via local opera- 
tions which do not change the entanglement of the state |7]. 

The main goal of our work is to focus on the idea of en- 
tanglement relativity Q, y[ as it can be applied in the field 
of gaussian states. As already pointed out, once a gaussian 
state is given in terms of its covariance matrix with respect to 
a given tensor product structure, we study how the entangle- 
ment properties of this state change when new tensor product 
structures are defined via number conserving modes transfor- 
mations. The latter will be defined in detail in the following 
sections. It is important to emphasize that this paper is not 
primarily aimed at studying gaussian states per se. The lat- 
ter are used, thanks to their structural simplicity, as an ideal 
arena to address and illustrate in a quantitative way the some- 
what abstract idea of entanglement relativity. It is important 
to stress the difference, with respect to the existing literature, 
of the conceptual status of the modes transformations that we 
are going to consider. Indeed, in most of the former studies 
(e.g., (alX modes transformations are physically enacted by 
optical devices (e.g., beam-splitters) that actually change the 
quantum state under examination. In this paper instead, the 
quantum state p of two bosonic modes, has to be thought as 
unchanged: it is just the way p is decomposed into subsys- 
tems i.e., modes, that is modified. The situation is somewhat 
reminiscent of the difference between "active" and "passive" 
rotations in elementary geometry; just as in that case, while 
the mathematical objects involved can be the same, the con- 
ceptual interpretation is quite different. 
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The paper is organized as follows. In Sect. (JIIJ the ba- 
sic ideas of entanglement relativity are reviewed. In Sect. 
(II I It theoretical background on gaussian states is provided. 
Sects (IIV A> and (II V B> are devoted to the analysis of pure 
and mixed gaussian states respectively. Finally Sect. (|V} con- 
tains some discussions and conclusions. 



elements are parametrized by finite-dimensional objects i.e., 
Nby N unitary matrices. 

III. THEORETICAL BACKGROUND 
A. Gaussian states 



II. ENTANGLEMENT RELATIVITY 

In the general setting |3] the local commuting algebras we 
are going to consider in the following are those generated by 
the mode creation/annihilation operators the global state- 
space considered here is the Fock space (total number of par- 
ticle unconstrained) that it is known to be isomorphic to the 
tensor product of harmonic oscillator state spaces. It is an el- 
ementary fact that, by using canonical transformations, one 
can move from a given set of modes to a new one defining 
a new tensor product structure i.e., new subsystems. One is 
then naturally led to analyze how the entanglement proper- 
ties of a given quantum state p depend on the choice of the 
set of modes and how they change when a different set of 
modes is considered. In particular, once p is given, one may 
want to know whether it exists a choice of modes i.e., a ten- 
sor product structure, with respect to which p is a separable 
state. In the negative case one may refer to such a state as 
absolutely entangled. Another natural task is to find the mini- 
mum and the maximum of the entanglement of p over the set 
of all possible tensor product structures (i.e, modes redefini- 
tions). Besides its theoretical interest this might have some 
practical relevance. For instance,the latter analysis allows one 
to identify for any given state the modes that correspond to the 
maximal entanglement and that can be fruitfully exploited for 
some quantum information processing protocol e.g., telepor- 
tation. 

More formally, if one denotes by T S TVS one of all 
the possible tensor product structures (TPS) and by Ef an 
entanglement measure e.g., negativity, associated to the TPS 
T one can define at least two natural "absolute" i.e., TPS- 
independent, entanglement measures 



E+(p) := sup E T (p), 

T^TVS 



E-{ P ) 



inf Er{p). 

TeTVS 



(1) 



The set of absolutely separable 1 18] and absolutely entangled 
states would then correspond respectively to {p / E + (p) = 0} 
and {p I E~(p) > 0}. As mentioned above a very natural 
question that arises concerns the relative weight of these two 
subsets in full state space. Of course, a more piratically mean- 
ingful definition of a TPS -independent entanglement measure 
would involve in Eq. Q a restricted set of operationally 
meaningful TPSs. Indeed, it should be clear that the defini- 
tions will in general involve optimization procedures over 
a vast set of possibilities. 

In order to make this program amenable to a quantitative 
treatment, in this paper we shall consider bi-partite gaussian 
states only. This allows us to focus on a subset of TVS whose 



Gaussian states are defined as those quantum states whose 
Wigner function is gaussian; they are completely defined by 
their second statistical moments. By grouping all the compo- 
nents Xi of the position and momentum operators in a single 
column vector such that Xi — Xi+i = pi, one can write 



-~(x*V-ix) 



(2) 



where W(x) is the Wigner function and V is the covariance 
matrix. The entries of the latter are exactly all the second 
order statistical moments of the components of the canonical 
operators (position and momentum). 

In general, in terms of the components Xi, the classical co- 
variance is defined as: < XiXj > — < X{ >< Xj >. Dealing 
with quantum operators the products appearing in the covari- 
ance expression must be taken in the Wigner-Moyal form, i.e. 
they must be symmetrized: < x.- L Xj >—< XiXj + XjX. L > /2. 
Since we work with gaussian states, the first moments (mean 
values) can be set to zero by means of displacement operators 
that do not change the entanglement properties iToll . There- 
fore V is built by simply setting: Vij =< XiXj >. In this 
way , the covariance matrix V has all the variances < xf > 
on the diagonal and the covariances < XiXj >, i ^ j off di- 
agonal. Due to the Wigner-Moyal symmetrization V turns out 
to be real and symmetric. In the following we will deal with 
bipartite systems: the identification of x±, X2 with qA,PA and 
of X3 , X4 with qB,PB, allows us to write V in a block form: 



V 



A C 
C* B 



(3) 



where A, B refers to subsystem A(lice) and B(ob) respec- 
tively. 

By using the covariance matrix formalism the Heisenberg 
uncertainty principle can be recast in the following compact 
way: 



V>-a, 



(4) 



where a is the so-called symplectic unity for two subsystems: 



1 CT 

-10 

1 
0-10, 



(5) 



Here, the uncertainty principle is expressed in its extended 
form; together with the variance relations, also the covari- 
ance relations are taken into account: Far-fa^] Var[x.,'] — 
Covix^Xj} 2 > 1/4. 
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Any covariance matrix V can be put in the so called stan- 
dard form V s by means of a suitable local symplectic trans- 
formation (see below). In V s the blocks have a diagonal form 



V s = 



a 


V o 



o 

a 


c_ 



\ 



(6) 



Since the transformation that gives V s is local, it does not 
change the entanglement properties of the state. 



B. Modes transformations 

A generic bosonic bipartite (gaussian) state can be de- 
scribed in terms of two harmonic oscillator modes. In a quan- 
tum system the operators accounting for the oscillation ampli- 
tudes are the annihilation/creation operators ai,a\. The two 
subsystems corresponding to the bipartition are the two oscil- 
lation 'directions'. We can imagine Alice and Bob sharing the 
system, Alice having one mode, Bob having the other one. 
According to what we have already said about the relativity of 
the entanglement, one state can be entangled if described by 
a particular pair of modes and can be separable if described 
by a different pair. If we now consider how the modes can 
be changed we see that, from a phase-space point of view, the 
modes transformations correspond to canonical transforma- 
tions of the conjugated (position and momentum) hamiltonian 
variables. The new variables must of course fulfill the same 
canonical commutation relations. The set of these canoni- 
cal transformations form the symplectic non-compact group: 
Sp(2n,TZ). 

A well known class of mode transformations are the 
so called Bogolioubov transformations: a — > b, 6j = 

J2j a ij a j + Pij a j> w i m hj — 1) —,n. For bipartite states 
n = 2. A special sub-class of Bogolioubov transformations is 
given by: 



(»,i = i,2) 



(7) 



where U is a unitary 2x2 complex matrix. This kind of mode 
transformations are those which conserve the number of par- 
ticles and are precisely the ones we are going to focus on in 
our analysis. These mappings form a group K(n) that coin- 
cides with the compact subgroup of Sp(2n, 1Z) of maximal 
dimension full . Each mapping @ induces a specific tensor 
product structure T S TVS. The goal of our paper is then to 
study how entanglement varies when the modes redefinitions 
belonging to K(n) are applied to a given state. In particular 
this means that in the remainder of the paper the ensemble 
TVS used to evaluate is in a one-to-one correspondence 
with elements of K(n). 

Each of the unitary and symplectic transformations Su E 
K(n), where U is defined in @, can be written as l2lll : 



S„ = 



X -Y 
Y X 



(8) 



where X and Y are respectively the real and the imaginary 
part of the unitary 2x2 operator U = X + iY acting on the 
modes fllll : 



X = 
Y = 



U + U* 
2 

U -U* 
2i ' 



The covariance matrix transforms according to l|8} as 

V = SuVSl. 



(9) 



(10) 



We end this subsection by discussing how one can represent 
a general element of K(n). In general a 2 x 2 unitary matrix 
U depends on four real parameters; thus a possible and useful 
parametrization of U is the following: 



U = 



sin tie 



cos 



0e ict> " 



cos Oe^ . , , 



The latter is the parametrization we use in our numerical cal- 
culations. 



C. Entanglement measure 

In continuous variables systems a proper measure of entan- 
glement that can be efficiently computed is logarithmic neg- 
ativity 15 fl^. Eoll . Following Ref. |6] logarithmic negativity 
can be obtained starting from the trace-norm of the partially 
transposed density matrix. This trace-norm can be computed 
from the so-called symplectic eigenvalues of the partial trans- 
pose of V with respect to the subsystem A, i.e. the eigenvalues 
of v~ 1 V Ta . The characteristic equation for ct~ 1 V Ta is: 

A 4 - (dct A + dct B - 2 dct C)A 2 + det V = (12) 

Thus the logarithic negativity can be computed as: 

^AT = i^max[- log 2 C2|A i |),0]. (13) 



IV. RESULTS 

In order to give a quantitative description of how the en- 
tanglement of a gaussian state changes under the action of all 
unitary transformations Sjj € K(n), we have performed dif- 
ferent kinds of numerical computations that now we briefly 
describe. 

/. Pure states We start our analysis by considering pure 
gaussian states whose covariance matrix is in the standard 
form V s ■ The modes transformations are generated by using 
a suitable discretization of the four dimensional parameters 
space defining U, see (II li . For each transformation (mode 
decomposition) we compute the logarithmic negativity of the 
corresponding covariance matrix (1101 . 

Obviously the covariance matrices in the standard form do 
not describe all the possible gaussian states. V s is useful be- 
cause it depends only on four parameters and has the same 
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entanglement of a large class of states. Given a generic gaus- 
sian state and its covariance matrix V it is easy to bring the 
latter in standard form by applying a suitable local symplec- 
tic transformation Hie Sp(2,TZ) © Sp(2,K) such that 
V s = LVL 1 . One can devise the following procedure to gen- 
erate a set of generic V^s starting from the V s : given a covari- 
ance matrix V s one can apply the full set of local symplectic 
transformations T £ Sp(2,TZ) ®Sp(2, 11); for each Tone ob- 
tains the matrix V = TVgT 1 that corresponds to a squeezed 
state with the same entanglement of V s and with a possibly 
different energy. Once a generic V has been generated, we 
apply the modes redefinition scheme SuVSjj and we evalu- 
ate the corresponding logarithmic negativity. The procedure 
is then repeated for many different V s and it is schematically 
displayed in Fig. Q. 







Fixed entanglement 


Sv 




/ 




V s 




/ 

Fixed energy 


T 






Fixed energy 


Su 



FIG. 1 : A visual display of the procedure used to generate and de- 
compose the gaussian states. 



We would like to underline that while the Sjjs transforma- 
tions (also referred in the literature as passive) are energy con- 
serving, the Ts are entanglement conserving and that in gen- 
eral SjjT TSu- It is thus important to stress that apply- 
ing the mode decomposition analysis to covariance matrices 
in non-standard form allows us to access a wider class of re- 
sults. In this respect the present work differs from what has 
been done in several studies on gaussian state entanglement, 
e.g. 1§L [T^ll . Those studies were indeed aimed at analyzing 
just entanglement properties of gaussian states, therefore the 
restriction to V s did not prevent the achievement of full gen- 
erality. Here our aim is to study the variation of entanglement 
for different bipartitions associated with unitary modes redef- 
initions, no matter if the initial covariance matrix V is V s or 
V. Thus the Ts transformations are merely used as a tool 
to define a more general set of new gaussian states; accord- 
ingly they do not have to be regarded as a new set of (number 
non-conserving) Bogolioubov transformations applied to the 
initial V s ." 

II. Mixed states. A great number of mixed states have been 
generated in a random way by drawing the parameters a, b and 
c± out of a uniform distribution over a well defined parame- 
ters range (see caption Fig. 0. Each (real symmetric) matrix 
is checked in order to verify whether it fulfills Heisenberg un- 
certainty principle; this allows one to select only the physical 
covariance matrices. The mode redefinition procedure is then 



applied and the maximal and minimal entanglement achiev- 
able is computed, see Eq. Q. These computations have in- 
volved covariance matrices both in standard and non-standard 
form. 

In the following sections we describe in details the results 
of our computations. 



A. Pure states 

The first set of gaussian states we take into accountwhose 
covariance matrix is in the standard form V s . In terms of V s 
they are described by only one independent parameter. In fact, 
with the usual notation (j§Jl, pure states correspond to covari- 
ance matrices with a = b (pure states are always symmetric), 
c± = ±V a 2 — 1, with the constraint a > 1/2. The pure gaus- 
sian states fulfill the Heisenberg uncertainty principle with the 
equal sign i.e., they minimize the uncertainty relations. We 
can consider two cases. In the first one c = 0: all covariance 
terms are zero and the corresponding states minimize the un- 
certainty relations in the usual form (AqiApi — 1/2). In the 
second case, c / 0, the pure gaussian states corresponding to 
V s are the squeezed states II 211 . It's not difficult to see that 
the non-squeezed states (c = 0) correspond to the isotropic 
oscillator case. 

The results of the numerical computations show that : 

i) there is always a unitary redefinition of the modes for 
which any pure state is separable i.e., E~ (p) = 0; 

ii) in particular, the non-squeezed states are absolutely 
separable: E + (p) = 0; 

Hi) in the case of squeezed states all unitary transformations 
lead to redefinitions of the modes (TPS) for which the 
value of the entanglement is always lower than value 
corresponding to the initial V 8 . 

A few comments are in order. The existence of a unitary 
redefinition of the modes for which the state is separable, see 
i), is not surprising. In fact this can be understood, by ele- 
mentary arguments, in different ways. One one hand one can 
resort to the situation in which two coupled oscillators can be 
decoupled by the classical normal modes decomposition. A 
suitable choice of the mode operators (either classical or quan- 
tum) always allows one to deal with two independent mode 
of oscillation. On the other hand, the state is described in 
the configuration space by the Wigner function W of Eq. @. 
If we consider W as a function of qA,qB and a plane par- 
allel to the qA,qB plane intersecting W(qA,qB), we get an 
ellipse whose axes lay on the bisectrices of the quadrants. In 
the isotropic case (\c±\ ~ 0), see ii), the ellipse reduces to a 
circle, therefore the modes redefinitions, correspondig to ro- 
tations in the plane, do not change the variances of qA,qB- 
When \c±\ > a covariance term appears and it is respon- 
sible for a non-vanishing excentricity of the ellipse, i.e. the 
gaussian Wigner function is 'squeezed' . A rotation (mode re- 
definition) of the reference axes always allows one to make 
them conincide with the ellipse's ones; therefore, it is always 
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possible to set to zero the covariance term, which is respon- 
sible for the entanglement between the modes. In our case a 
7r/4 rotation reduces the entanglement to zero (see Fig. In 
the q A (B)PA(B) Plane instead, since Cov[q A{B) ,p A(B) ] = ( 
i.e., the blocks A and B in Eq. Q are diagonal), the argument 
of Eq. simply identifies a circle. 

We now turn to examine point Hi), In Fig. |2]the logarithmic 
negativity Jl 31 of different pure sates p is plotted as a function 
of 6 (cj), ft, (j)" in Eq. lO are kept fixed). 



1 1 1 1 1 1 1 
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FIG. 2: Logarithmic negativity as function of the parameter 8 for 
different pure states. The picture clearly shows a symmetric periodic 
behaviour. The point of zero entanglement is reached for 8 = n/4 
while the maxima are in correspondence of 8 — 0, 8 — tt/2. 

The maximal value of EV i-e-, E + (p), is achieved in corre- 
spondence of 9 = (modes swapping) and 9 = tt/2 (identity 
transformation). Any other value of 9 produces a lower value 
of Ej\f. In particular, for every state p, there is a value of 9 
for which E_\f = 0. The results, that would not qualitatively 
change with a different choice of <j), <j>' , (/>", can be understood 
as follows. In case of symmetric states the maximal entangle- 
ment is achieved in correspondence of V s because, in this case, 
all the quantum correlations between the operators are set to 
zero except those involving same-type operators between dif- 
ferent subsystems. Every unitary modes redefinition gives rise 
to correlations between conjugated operators belonging to the 
same subsystem, while it decreases the correlation terms and 
consequently the entanglement between the two subsystems . 
Similar results were obtained in the context of optical passive 
transformations on gaussian states by Wolf et al.|8] and they 
are in agreement with the fact that the state described by V s 
corresponds, at given energy, to the maximum entanglement 

Q. 

We now turn to discuss the states whose covariance matrix 
V = TV s T l is not in standard form. It is important to realize 
that, at variance with the former studies of gaussian state en- 
tanglement, in the present setting the passage to the standard 
form is not harmless. In fact in order to cast a generic covari- 
ance matrix V into its standard form V s , one has to perform 
local operations L ( e.g., local squeezing) that are outside the 
class of the modes redefinitions here allowed l|8}. Since these 
latter do not in general commute with the local operations L 
required to get V s one obtains, as we will now see, rather dif- 
ferent results depending whether the standard form is consid- 



ered or not. 

The numerical computations, see Fig.|5J show that: 

iv) as described for the previous case, it is always possible 
to find a bipartition in which the state is separable i.e., 

E-(p) = 0; 

v) it always exists a redefinition of the modes that can in- 
crease £V- 

Fig. 10} shows the behaviour of the logarithmic negativity as a 
function of one of the local symplectic parameter (the squeez- 
ing operation) -q : at — > diag(rja, r/~ 1 a, rja, ?7 _1 a) and of one 
of the global unitary transformation parameter 9. 




e ti/2 



FIG. 3: The main difference between the behaviour of the pure state 
with V s (dotted line) and the one with V (solid line) is the possibility 
for the latter to increase its logarithmic negativity 




FIG. 4: A 3-dimensional plot of the logarithmic negativity as func- 
tion of sind and r\. It is possible to note a line laying on the ground 
plane (dark valley) testifying the zero-entanglement modes redefini- 
tion for each initial state. Also it is worth noting the increasing value 
of as r\ gets different from 1, for some [7^1. 
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B. Mixed states 

The numerical analysis of mixed states leads to other inter- 
esting results. As in the previous section we start by con- 
sidering the the covariance matrices in standard form. We 
first look at the symmetric (a — b) subset of mixed states 
depending only on three parameters. If we set, as an addi- 
tional condition, c + = — c_, we have states that behave in a 
way similar to the pure ones. In fact Ejy is maximal in corre- 
spondence of the identity transformation and it decreases for 
all entangling [7^1. The difference consists in the fact that, 
instead of a single point, there is a whole interval of values of 
9 for which the logarithmic negativity vanishes, as displayed 
in Fig. |5] The further one moves away from pure states (ps) 
(\c±\ < |c± s | = vV-1) the wider is the plateau. When the 
plateau covers all the range of variation of 9 the states become 
absolutely separable. The case \c±\ > \c± \ would be the only 
one for which there is no point of zero entanglement, but it has 
to be discarded since the corresponding matrices do not fulfill 
the Heisenberg principle. For simmetric mixed states such 




FIG. 5: The plot of the logarithmic negativity as a function of 8 for 
three different symmetric mixed states shows a qualitatively similar 
behaviour for the three states. The width of the vanishing entangle- 
ment plateau is related to difference between \c± | and |cj | 

that V 7^ V s the entanglement calculations give results that 
are similar to the ones obtained for the pure states case (see 
Fig|6): unitary redefinitions of the modes can increase the en- 
tanglement. In addition, as described above, E^ can vanish 
in correspondence of a whole interval of mode redefinitions. 

Let us now turn to analyze the case of non-symmetric mixed 
states, starting as usual from those whose covariance matri- 
ces are in standard form. A great number of mixed states, 
in the defined range, turn out to be absolutely separable i.e., 
E + (p) = 0. FigQ\ shows the amount of absolutely separable 
states with respect to the amount of randomly generated phys- 
ical states. Although the chosen range is not representative of 
the whole mixed gaussian states space i.e., different choices 
would have led to different results, it allows us to make some 
general remarks. 

The high separability rate can be understood by using the 
Simon criterion 1 9] : it is sufficient for c± to have the same sign 
to achieve separability. This set of absolutely separable mixed 
states includes all the classical states defined in Ref. fl4llT3l . 




e jt/2 



FIG. 6: Logarithmic negativity as a function of 9 for a mixed sym- 
metric state (solid line) with covariance matrix in non-standard form, 
compared to the one of the corresponding state with V 3 (dotted line). 



On the other hand, when the mixed state are not absolutely 
separable there always exists, as in the pure states case, a non- 
void set of unitary modes redefinitions for which the entan- 
glement is zero. This result, together with the previously ob- 
tained for pure states, allows one to claim that no absolutely 
entangled bipartite gaussian states exist. Every gaussian state 
can be regarded as separable, it is just a matter of choosing 
the right operators/modes redefinition. The maximum entan- 
glement is reached by pure states, while mixed ones have, in 
general a low logarithmic negativity. 




FIG. 7: Statistical results for gaussian mixed states with covariance 
matrix in standard form (A) and V ^ V s (B). The parameters range 
is the following: Vij G [-1, 1] for i j while V» G [0.5, 1.5]. 
The rate of absolutely separable in first case (A) states is about 94% 
over 7746 generated states, while in the second case (B) the rate of 
absolutely separable states is lower (about 42% over 4692 states). 

We finally consider the case of non-symmetric mixed states 
for which V ^ V s . The new statistics (FigQi), within the 
same parameters range, shows that, in comparison with the 
V = V S case, the set of states that have nonzero entanglement 
for some bipartition is far larger. This is due, as in the pure 
states case, to the fact that entanglement can be increased with 
a proper choice of mode redefinition and it is true even for 
initially separable states. We remind that a different choice of 
the range would lead in general to different percentage rates. 
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V. CONCLUSIONS 

The main purpose of this paper is to show how the con- 
ceptual issue of entanglement relativity i.e., the crucial depen- 
dence of entanglement on the choice of subsystems, can be 
addressed in a quantitative fashion for a specific, yet quite im- 
portant, class of quantum states. We have studied how the 
entanglement of a two-modes bosonic gaussian state changes 
when the modes are redefined through number-conserving 
unitary transformations. This redefinitions have to be seen 
analogously to "passive" rotations in elementary geometry: 
the quantum state is unaffected, it is just the way it is split 
into subsystems i.e., modes, that changes. One obtains a set 
of bi-partitions parametrized by 2 by 2 unitary matrices; this 
makes the choice of the set of possible bipartitions, typically a 
huge and rather ill-defined one, easy to be dealt with analytical 
and numerical techniques. 

We have found that with a suitable unitary redefinition of 
the modes every gaussian state can be regarded as separable. 
This means that there are no absolutely entangled bipartite 
gaussian states i.e., for the chosen class of tensor product 
structures the function E~ (p) (see Eq. Q is identically van- 



ishing. Moreover a large class of mixed gaussian states can 
be considered as absolutely separable, E + (p) = 0, meaning 
that no unitary redefinition of the modes can 'entangle' these 
states. 

This kind of entanglement variation analysis can be ex- 
tended to a larger set of Bogolioubov transformations: those 
that do not conserve the particle number and that are able to 
induce further bipartitions 1 22] . 

Of course even upon this extension the main result about 
the non-existence of absolutely entangled gaussian states will 
remain valid, whereas the set of absolutely separable states 
might become smaller. Finally, recent papers Il6l [nil have 
addressed the problem of dealing with gaussian states even in 
fermionic systems. This naturally suggests the possibility of 
extending our investigations to the fermionic realm. 
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